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We develop the theory for the Adiabatic Raman Photoassociation (ARPA) of 
ultracold atoms to form ultracold molecules in the presence of scattering resonances. 
Based on a computational method in which we replace the continuum with a discrete 
set of "effective modes" , we show that the existence of resonances greatly aids in the 
formation of deeply bound molecular states. We illustrate our general theory by 
computationally studying the formation of ^^Rb2 molecules from pairs of colliding 
ultracold ^^Rb atoms. The single-event transfer yield is shown to have a near- 
unity value for wide resonances, while the ensemble- averaged transfer yield is shown 
to be higher for narrow resonances. The ARPA yields are compared with that 
of (the experimentally measured) "Feshbach molecule" magneto-association. Our 
findings suggest that an experimental investigation of ARPA at sub-/iK temperatures 
is warranted. 



I. INTRODUCTION 



Adiabatic Raman Photoassociation (ARPA) of ultracold atoms was introduced [IHS] as 
a practical way of producing ultracold diatomic molecules in their ground electronic and 
vib-rotational states. As illustrated in Figjl} the method consists of photoassociating two 
colliding atoms by two ("dump" and "pump") laser pulses that are mutually coherent and 
partially overlapping in space and time. As in three-level "Stimulated Raman Adiabatic 
Passage" (STIRAP) |1HS], one uses the "counter-intuitive" pulse ordering [1] in which the 
"dump" pulse, connecting the final bound state to an intermediate excited bound state, 
precedes the "pump" pulse, connecting the continuum to the latter state. In this way one 
executes a smooth "Adiabatic Passage" from an (ultracold atom-scattering) continuum to 
deeply bound molecular states 

The main problem with the above approach is the small "Franck- Condon" overlap factors 
between the intermediate bound state and the continuum. The introduction of a ("Fesh- 
bach") resonance which is expected to better overlap with the intermediate bound state 
can alleviate this problem [HI HO]- As we show below, a Feshbach resonance also induces 
an important dynamic effect of prolonging the lifetime of the spatial region (the "Franck- 
Condon window") in which photoassociation occurs. In this way a larger fraction of the 
atom pairs becomes available for photoassociation. This aspect has not been considered in 
previous works on photoassociation via Feshbach resonances [TU], where it was concluded 
that wide resonances are more effective in promoting photoassociation than narrow ones. 
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FIG. 1: A schematic display of the Adiabatic Raman Photoassociation (ARPA) process. Left 
panel: Atoms colhding in the near-threshold energy range are excited by the "pump" laser to 
the vibrational state |2) on an excited electronic potential. The latter is coupled by the "dump" 
laser to the deeply bound target state Right panel: The same as on the left panel, with the 
continuum-bound couplings modified due to presence of a Feshbach resonance. 

While we agree with Ref. |ilOj concerning the outcome of a single pair collision, we differ 
in our conclusions regarding the ensemble averaged yields: Based on our calculation of the 
number of recombining atoms per laser pulse-pair, we find that narrow resonances are in 
fact more effective than wide ones, because by prolonging the lifetime of the Franck Condon 
window, narrow resonances allow for more recombination events to occur. This mechanism 
more than compensates for the smaller energetic widths of the narrow resonances. 

The structure of this paper is as follows: in section II we develop the working equations 
for the ARPA process, based on representing the continuum as a discrete set of "effective 
modes" [H]. In section III we present calculations of the ARPA population dynamics for a 
single pair of ultracold ^^Rb atoms, and demonstrate how the FC window lifetime and spatial 
extent are being extended by the resonance. In section IV we compare the yield of ARPA 
with that of magneto-association [T2HT6] in which a pair of atoms in a Feshbach resonance 
are transformed into a stable molecule by sweeping over an external magnetic field, thereby 
pushing the resonance energy to lie below the molecular dissociation limit. We show that 
the two schemes lead to similar scaling of molecular production yield, but that ARPA is 
expected to be more efficient. Finally in the Appendix we show that the action of the pulses 
is tantamount to a quantum projective measurements on the initial continuum wave packet 
[21 E]- Concluding remarks are provided in section V. 
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II. THEORY 

A. The basic formulation 

As illustrated in FigjT| ARPA involves a A-type level structure, similar to 3-state STIRAP, 
in which two, mutually coherent, temporally and spatially overlapping, laser pulses induce 
adiabatic passage from a molecular continuum (representing two colliding atoms) to the 
target bound level using an excited bound state |2) as an intermediate. The Hamiltonian 
of the system is written as (in atomic units), 

H = Hq — 2fi e„(t) cosWnt (1) 

n=l,2 

POO 

where i^o = + ^2|2)(2| + / E\E){E\dE, (2) 

is the "material" Hamiltonian. The bandwidth of the pulses ranges from being of order of 
100/iK, down to tens of nK, which, when compared to the vibrational energy separation, 
makes it valid to include no other bound states than |1) and |2). The second term in Eq.([T]) 
describes the interaction of fi, the transition dipole moment, with the "dump" [n = 1) and 
"pump" {n = 2) laser pulses, whose respective amplitudes and central frequencies are e„(t) 
and Un- We tune U2, the pump center frequency, to be in near resonance with U2,e, the 
transition frequency between the intermediate state |2) and the continuum states \E); and 
uJi, the dump center frequency, to be in resonance with W2,i, the intermediate-to- final-state 
transition frequency. We assume that the laser fields do not vary significantly over the range 
of atom-atom distances in which photoassociation occurs, thereby justifying the elimination 
of the spatial variation of the fields. 

As state |1) we choose a deeply bound rovibrational level belonging to the ground elec- 
tronic potential. Given this choice, the intermediate state |2), taken to belong to an excited 
electronic potential, is chosen to be a vibrational state that overlaps well with the |1) state. 
The main feature of the continuum we explore here is the embedding of a (Feshbach-type) 
resonance, leading to a sharp energy dependence of continuum-bound transition-dipole ma- 
trix elements /zs.e = W\E) 

Expanding the time- dependent system wave function in the material basis set, 

POO 

I^W) = Yl biit)e-'''''\t) + / dEbEit)e-'^'\E), (3) 

1=1,2 ''^th 

we obtain, using the time-dependent Schrodinger's equation i^|\Ef(t)) = i^|\E'(t)), the or- 
thonormality of the material states and the Rotating Wave Approximation (RWA), that, 

hit) = tniitMt) (4) 

POO 

hit) = ^n^it)b^it)-Tfb2it) + l / nEit)bEit)e'^^'dE (5) 

J Efh 



bEit) = iniit)b2it)e-'^^\ (6) 
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where Eth is the continuum threshold energy, Ae = E2 — E — L02 and Ai = E2 — Ei — ui 
are the detunings of the pulses, and T/ is the spontaneous decay rate of state |2). There are 
two Rabi frequencies in the problem, Qi(t) = ei(t)/i2,ie*^i*, and QE(t) = ^2(t)fi'2,E- 

Eq. ^ for bE^t), representing a continuous set of equations for the continuously varying 
E, are numerically difficult to solve. We therefore eliminate Eq. ^ by integrating bE{t) in 
time. 



bE{t) = bE{0)+t f dt'QUt')b2{t')e-'''^'' , 
Jo 

and substitute this formal solution into Eq. ^ to obtain that, 

POO 

b2{t) = ini{t)bi{t) - Tfb2{t) + I / nE{t)bE{Q)e'^'''dE 



(7) 



/oo r /"t 
/ e2{t')b2{t')e~^''-'' dt' 
■00 . Jo 



Eth 

iAst 



dE. (8) 



By defining fsource{t), the source function, and F{t—t'), the spectral auto-correlation function 

as 00 

fsourceit) = / fis(t)6i,(0)e*^^*rfE, (9) 
J Eth 

POO 

F{t -t')= / \ii2,E?e'^^^'-''^dE, (10) 

•J Eth 

we can transform the above (continuous) set of differential equations to a set of two integro- 
differential equations, 

61 (t) = l^ll{t)b2{t) 

b2{t) = I^l{t)b^{t)-Vfb2{t)+ ^fsourceit) -e2{t) [ e2{t')b2{t')F{t-t')dt'. (11) 

Jo 

The threshold energy Eth in the source function will later in our analysis be taken as — cxd, 
since the function fe_B(0) is zero near the collision threshold of the ground electronic potential, 
refiecting the density of states at zero kinetic energy. 



B. The effective modes expansion 



The simplest solution of Eqs. (11) is obtained by the "fiat continuum" or "slowly varying 



continuum approximation" (SVGA), according to which, whenever ^2,e varies sufficiently 
slowly with energy E we replace it by its value at some average energy E. In this way the 
spectral auto-correlation function of Eq. (10) is reduced to F{t — t') = 2n\^2E\'^^{'t' ~ t)y 



and the integration in Eq. (11) is eliminated. Given this approximation, the dynamical 



equations assume, in matrix notation, the form, 

d 



dt 



b = ■ b + it, 



(12) 
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where b(t) = {bi{t),b2{t))^ , isourceit) = (0, fsourceit))'^ , with ^ designating the transpose 
operation. The Hamiltonian matrix is defined as, 

H = I ^re//(t)) ^'^^ ^'"^^^ ^ ^\l^2,E\'4it). (13) 

A detailed discussion of the solutions under SVGA was made in Refs. [IH3]. 

The SVGA is however invalid when collisional resonances are embedded in the continuum, 
because in that case fi2,E changes rapidly near the resonance energy [9l [17] . In order to treat 
this case we parametrize fi2,E as [TT] 

s=l ' 

where /x^ represents the electronic transition dipole moment; - the FuU-Width-at-Half- 
Maximum (FWHM); and Es - the centre-position of each s resonance. This form is capable 
of approximating well both wide and narrow resonances [101 US HE] • As will be seen below, 
the above parametrization allows us to greatly simplify both the analytical theory as well 
as the numerical propagation of the dynamical equations. 

With the expansion (14), the auto-correlation function F{t — t') in Eq. (10) can be 
evaluated analytically as, 

M 

F{t-t') = Y,aJt{t)f;{t') (15) 

s=l 

with 



Using Eq. (15) we now define [TT] the "effective modes" variables as, 

B-{t) = i [\2{t'Mt')f;{t')dt', (17) 

Jo 

using which, we transform Eqs. (9) and (10) into, 

bi{t) = ini{t)h2{t) (18) 

M 

h2{t) = iQi{t)bi{t) - TMt) + ifsource{t) + ze2 (t) J] a (t)5; (t) (19) 

s=l 

B-{t) = ie2{t)f;{t)b2{t), s = l,...,M. (20) 



In this way the original set of continuous equations for &£;(t) is replaced by a discrete set 
of equations for (t) . We can further simplify the structure of the equations by defining 
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Bs{t) = a/ as/27rf^{t)Bg (t), and mode-dependent Rabi frequencies, = ^2{t)\/27m 



With these definitions Eqs. (20) assume the form 



hit) = tni{t)b2{t) 
62 (t) = ini{t)bi{t) 



^Mt) + ifsourceit) + ^ (t)5, (t) 



B,{t) = -tXsBs{t) + ini'\t)b2{t). 

Writing these equations in matrix notation, we have that, 

d 



dt 



b = zH • b + ifg, 



(21) 
(22) 

(23) 



(24) 



where 



b(t) 



Bi{t) 

B2{t) 

V '■■ J 



.it) 



\ 



\ 

source (t) 




■■ J 



f ni 



and H 










-Xi 








-X2 



(25) 
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The "effective modes" amphtudes Bs{t) thus appear equivalent to some extra bound states 
of energies Eg that are coupled by the Rabi frequencies ^2\t) to state |2), with detuning 
ES—E2+0J2 and decay rates Ts/2 as contained in Xs- The non-Hermiticity of the Hamiltonian 
is due not just to the decay of the effective modes, appearing as the imaginary part of Xs, but 

(s) 

also to the Rabi frequencies fig 5 which are in general complex numbers, since the definition 
of as involves a summation over s', namely the effective interaction between overlapping 
resonances. 

Eq. (24) resemble (multi-state) STIRAP [19] with fli{t) and fi2*^(t) coupling respectively 



|1) with 1 2), and |2) with each of the s effective modes (Fig. |2]). We note however that the 
transfer dynamics differs in a significant way from ordinary STIRAP in that here the initial 
population does not reside in the effective modes, which get gradually populated. We can 
see this most explicitly for a single resonance, for which the dynamical equations assume 
the form, 



61 (t) = ini{t)b2{t) 

hit) = iniit)bi{t) ~Tfb2{t) + tni'\t) 

Bi{t) = ~ixiBi{t)+t9!i\t)b2{t). 
By re-defining B{t) = fsouTce{t) /^2^\t) + Bi{t) we obtain 

bi{t) = iniit)b2{t) 

hit) = ini{t)bi{t) - TMt) + iQ^2 \t)B{t) 



fsource{t)/n^2\t)+B,it) 



B{t) = -tXiB{t)+tQ^^\t)b2{t) + 



^Xlfsource{t)/n^2\t) + f'sourceit) 



(26) 
(27) 

(28) 

(29) 
(30) 

(31) 
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FIG. 2: Left panel: Three bound states STIRAP. The population which starts in the initial 
bound state |3) is transferred into state |1) by following the evolution of the "dark" field-dressed 
state. The process avoids population loss due to spontaneous emission from state |2) because the 
latter remains unpopulated. Middle panel: ARPA via a collisional resonance. The population 
gradually feeds the resonances as the continuum wave packet (shaded area) arrives (at time io)- 
Right panel: The resonance-dominated continuum of the middle panel is replaced by a single 
(several) effective mode(s) with decaying amplitude(s) Bi(t) {Bs{t)). 

where f'sourcei'^) = 72!^ f^^'^Ef^2,EbE{0)^^^'^^dE. We see that here it is the terms in the 
square bracket that populates the effective mode (see Fig. |2]). 

Another contrast with three state STIRAP is the possibility of leakage of population 
from the "dark" state. In three states counter-intuitive pulse ordering adiabatic passage 
[1], the population resides initially in the adiabatic "dark" state, which is a superposition 
of the initial and target states only. In that case, the adiabaticity of the pulses guarantees 
the completeness of the transfer from the initial to the target state, leaving the intermediate 
state unpopulated at all times. Because in our case the effective mode gets populated 
in a gradual fashion, the system wave function may contain non-negligible contributions 
from other ("bright") adiabatic states. These "bright" states have a small overlap with the 
intermediate state |2), causing population to be lost via spontaneous emission. 



III. COMPUTATIONS OF RESONANT PHOTOASSOCIATION 

We view the entire ARPA process as a statistical average over collisions between indi- 
vidual pairs of colliding atoms. As explained above, we represent such pairs by a set of 
(spatially extended) coherent wave packets arriving at the Franck Condon region at differ- 
ent times. The choice of coherent wave packets (rather than plane waves) as the basis of 
our our computations is merely a result of our wish to work with normalizeable states. 

In this section we examine the above formulation by performing a set of computations 
on the resonantly enhanced photoassociation of ultracold ^^Rb atoms to form ^^Rb2 in its 
ground vibrational state. In keeping with our view of the process we divide the computations 
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Energy [^K] 

FIG. 3: A single resonance (in the effective mode expansion) is used for the computation, with 
magnitude |/i2,£;| = \^J'res^res/ — Ej-es) + ^r^es]! as a function of energy, shown in this figure 
for various choices of T^es- The resonance height has numeric value fj-res = 300a. u. (depending on 
the actual physical system, an even larger height value can be used, which favours a lower cost of 
laser amplitude required). The peaking shape of the resonance represents the enhancement of the 
FC factor. Centre of the resonance is Eres = lOO^K. 

into two parts: (A) population transfer for each colliding pair, and (B) population transfer 
of the thermal ensemble of pairs of colliding atoms. 



A. The single collision transfer yields 

Following the model of Refs. [21 13], we consider a pair of '^^Rb atoms colliding on the 
ground electronic potential. We assume that at t = 0, chosen to occur before the onset of 
the pulses, all the population resides in the continuum wave packet and none in states |1) 
or |2). The shape of |^(0)) = J^^^dEbE{0)e~''^^^\E), the initial continuum wave packet of 
Eq. (|3|, with Eth being the lower energetic limit (which is extended to — oo), is determined 
by the &£;(0) function, chosen here to be an energetically-narrow Gaussian [IHS] 



{E - Eof 



251 



i{E - Eo)tu 



(32) 



with ^0 = 70/iK and Eq = lOO/zK. With this choice of parameters, the wave packet temporal 
peak occurs at to = 1-2/is [HE]. The pair of pulses spectral widths are then chosen to 
have a good overlap with the energetic spread of the atomic wave packet. This requirement 
translates in the time domain to pulse durations. 

The scattering continuum is assumed to contain a narrow resonance, whose shape is given 



by Eq. (14). Figure |3j shows the shape of a resonance centred at Eres = Eq = 100/iK for 3 
different widths. Since the (Feshbach) resonance position, can be tuned (magnetically), we 
can optimize the transfer by tuning Eres to be always equal to Eq. In this way one achieves 
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the maximal FC enhancement (as confirmed by the numerical calculation presented in Figure 
5] bottom panel). 

The bound-to-bound matrix element is chosen to have a numeric value /j.21 = O.OOSla.u. 
This numeric value can be different depending on the actual experimental set-up, but our 
results only depend on the Rabi frequency fii(t), which is proportional to the product of this 
bound-to-bound matrix element with the dump pulse amplitude. A increase (decrease) of 
the matrix element translates into a proportional decrease (increase) in the laser amplitude 
required. The spontaneous decay rate from level |2) is Fj = (30ns)~^ The central 

frequency of the dump pulse is chosen to coincide with E2 — Ei, and the central frequency 
of the pump pulse - to coincide with E2 — Eq. The field amplitudes ei(t) and 62 (t) are taken 
as Gaussian functions, peaking, respectively, at 1.05/is and 1.55/xs. The duration of both 
fields is 0.22/is, and their peak intensity is 3 x 10^ W/cm^ [2J. 



With these specifications, Eqs. (21 23) simplify to yield 



6i(t) = tni{t)b2it) (33) 

hit) = ^^]l(^)6l(^) - vMt) + z/™(t) + ini'\t)B,{t) (34) 

Bi{t) = -{Tres/2)B,{t) + ini'\t)b2{t) (35) 

where the Rabi frequencies are ^i{t) = ei(t)/i2,i and 0.^2 ^ = e2(t)/ires\/27rTV^/2. Notice 
here that the FC factors contained in /i2,i and ^res always appear as product with the field 
amplitudes ei^2(t). So the intensities of the dump and pump fields really are determined 
by the respective FC factors between the intermediate and target states, and between the 
continuum and the intermediate state. An enhancement on either of the the FC factors 
will result in the same order decrease in the laser amplitude (square-root of the intensity) 
needed. 

In Fig. |4]we display the results of numerically integrating the equations for 61 (t), 62 (i) 
and Bi(t), given that &i(0) = 62(0) = -Bi(O) = 0. We plot the populations of states |1) 
and |2) for a wide (lOO/xK) resonance and a narrow (6/iK) resonance. The most striking 
feature of this plot is that the wide resonance gives rise to an essentially complete population 
transfer (> 90%), while the transfer probability via the narrow resonance is only ~ 23%. 
Since the target state is the ground state, no loss of population can occur after a single 
photoassociation event. Loss of population is however possible when subsequent collisions 
with the gas of atoms and/or the action of subsequent pulses are considered. As shown 
in the lower panel, due to the adiabatic nature of the process and the "counter-intuitive" 
pulse ordering, the population of the intermediate level |2), remains very low, even while 
the pulses are on. 
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FIG. 4: The target state and intermediate state populations as a function of time for two resonance 
widths. Top panel: the pump and dump field amplitudes. Middle panel: The target state 
probability The transfer yield is 90% for the wide (100/iK) resonance, but only 23% for 

the narrow (6/iK) resonance. Bottom panel: The intermediate state population |62(i)P- (Notice 
the large difference in the vertical scale relative to the middle panel.) 



B. Ensemble transfer yield 



In agreement with Ref. [TO], we have shown in sub-section A that for each event the 
transfer yield via a wide resonances is greater than that of a narrow resonances (Fig. 5 
top panel). The situation is however different for an ensemble of colliding atoms, where, as 
we show below, the transfer yield of the narrow resonances is greater. The reason is that 
for narrow resonances the number of colliding pairs which can react to the light is greater, 
essentially because for narrower resonances one can work with narrower pulse bandwidths, 
hence longer pulses. The increase in the number of effective collisions occurring during the 
increased pulse durations, more than compensates for the reduction in the individual event 
transfer yield. 

An alternative way of viewing this effect is to examine the role of the discrete effective 
modes which replace the continuum in our theory. These modes are to all intents and 
purposes resonances [20] . The only difference between the modes and scattering resonances 
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FIG. 5: Top panel: The target population, \bi{t — )• oo)p, for different resonance widths. Bottom 
panel: \bi{t — t- oo)p as a function of the centre of the resonance Eres, for £"0 = 100/iK; this shows 
transfer is optimal when the centre of the resonance coincides with the central energy of 6e(0). 



is that the effective modes do not originate from a real bound state embedded in a continuum. 



Thus, as clearly seen in Eq. (35), the rate of de-populating an effective mode is proportional 
Tres - the resonance-width of that mode. Hence narrower resonances, corresponding to 
smaller rates of depopulation, increase the interaction times of the effective modes with the 
intermediate level |2), thereby prolonging the duration of the Franck-Condon window. 

In Fig. |6]we examine these trends in a quantitative way by displaying f^(t), the field 
normalized source term, given as. 



f^it) = f,ource{t)/t2{t) 



lX2,EbE{0)e'^'''dE 



(36) 



for resonances of changing widths. Clearly in evidence is the prolonged duration of f^{t) 
when switching to narrower resonances. 

The temporally stretched population source is also beneficial when we consider the ac- 
tion of a pulse pair that is delayed relative to Iq, the arrival time of the incoming wave 
packet. Figure [7] shows the transfer efficiency as a function of such delay times for 3 dif- 
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FIG. 7: The \bi{t — )• oo)p transfer yield for different resonance widths as a function of the 

6t = to — tp delay time, where tp = 1.2/is is the pulses' overlap peak time, and Iq is the incoming 

wave packet peak time. 

ferent resonance widths. For a narrow resonance, despite the drop in the peak value, the 
single collision transfer efficiency remains large for longer times. This means that atom pairs 
which started their collision at an earlier time can still be transformed into bound molecules 
with non-negligible probability. 

In order to obtain the delay times averaged molecular production yield for an atomic 
ensemble we need to calculate the area under the transfer- yield curves of Fig. [7} Figure [8] 
displays the dependence of the delay times averaged yield for various resonance widths. We 
first note that the yield changes relatively slowly for resonance width larger than lOOOyuK. 
This is because in this case the resonance width by far exceeds 6q, the energetic spread of 
the initial atomic ensemble and we approach the fiat-continuum limit. As the width of the 
resonance drops to a few /iK, the molecular production yield rises to a maximum value, but 
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FIG. 8: The time-averaged molecular production yield for a single pulse pair for an atom ensemble 
at 100/xK, calculated by integrating the delay plots of Fig. [7] 

drops significantly due to spontaneous emission for yet narrower resonances. Thus there 
exists an optimal resonance width for which the molecular production yield is maximal. 
Comparing the optimal molecular production yield, obtained for a (narrow) resonance value 
of ~ 8/iK, with the yield in the flat-continuum limit, we see an improvement factor of 1.56. 



C. Scaling behaviour with ensemble temperature 

We now explore, as was done in Ref. [1], how the process varies as the average ensemble 
energy, Eq, and energy spread, Sq, are scaled down by a factor of s > 1, i.e., Eq -> and 
5o — ^- In ref. [1], we showed that the equations are invariant to this scaling provided the 
peak time was scaled up by the same factor to sto, and the initial wave packet amplitude 
is scaled as bE(0;E,to) — )• y/sbE {O; ^,sto). We now consider the effect, in addition to 
the above, of scaling the resonance shape as fi2,E{E,Tres) 1^2,e (f^^^)- order to 
match the spectral profile to the scaled 6£;(0), we need to scale up the centre frequencies 
and durations of the two pulses by the same s factor. Since we can choose the intensity 
(amphtude) of the fields, we scale ei(t) — )■ ^^^^ and 62 (t) — )■ [1] As a result, the source 
function scales like /source ^O; Tres) ■ fsource{st, sto, 1^). The above scaling leaves 



the dynamical equations (Eqs. (33 35)), essentially unchanged, except for the spontaneous 
decay rate which cannot be scaled. As we scale the relevant times by a factor of s, the 
deleterious effect of the spontaneous emission becomes more and more pronounced. 

In Fig. [9] top panel we display the dependence of the single collision photoassociation 
yield on the spontaneous decay rate. Note that a change in one order of magnitude for the 
spontaneous decay rate only affects our results negligibly. When we compare the results 
to those displayed in Fig. [5] top panel, where the transfer yield is plotted as a function 
of the resonance width, we see that the transfer yield is not greatly affected at ensemble 
temperatures of a few fiK to a few lO^/xK. As the ensemble temperature goes down by three 
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FIG. 9: Top panel: The single collision photoassociation yield \bi{t — )• oo)p as a function of 
resonance width, for four different values of the spontaneous decay rate. Bottom panel: The 
same plot at stronger laser intensities. The photoassociation yield now becomes insensitive to the 
spontaneous decay rates. 

orders of magnitude, the single collision transfer efficiency goes down too, by ~ 67% of its 
original value. The effect is more pronounced for narrower resonances, because the longer 
interaction times enhance the effect of the spontaneous decay. However, as shown in the 
lower panel of Fig. |9} it is possible to combat the effect of spontaneous decay at very low 
temperatures, e.g. in nK range, by increasing the amplitude (intensity) of both laser fields. 

One can summarize these results by saying that the optimal resonance width is always 
~ 8% of the ensemble temperature, and that the optimal molecular production efficiency by 
ARPA is ~ 56% higher than that of the wide resonance (flat continuum) case. 



D. Thermalization 

At long times the cumulative action of many pulses can change a region (or regions) 
in phase space corresponding to the recombining atoms, thereby affecting the initial wave 
packet amplitude 6£;(0). However, an atomic ensemble can thermalize sufficiently fast, on 
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the order of milliseconds, to yield the typical atomic trap setting of lOO/xK temperature and 
10^^ /cm^ density [26]. This means (depending on the repetition rate of the pulses), that 
after a few thousand fis pulse-pairs, the atomic ensemble can thermalize back to its original 
phase-space distribution, re- validating the ensemble- averaged form we used for 6£;(0). 

According to previous estimates [H |2], the total number of pulse pairs needed to transfer 
an entire atomic ensemble of density 10^^ /cm^ is around of 10''. Therefore a few thousand 
pulses is indeed a very small fraction of total number of pulses needed, and the thermalization 
is fast comparing to the ensemble size molecular conversion time. 

One is also inclined to pose the practical question of how to hide the newly formed 
molecule in state |1) from subsequent pulse pairs. In accordance with more detailed discus- 
sions in Ref PQ 12], this can be done, for example, by allowing the newly formed molecules in 
state |1) to "leak" away from the laser focus, which is possible because they react differently 
from the atoms to the confining laser frequency. A molecular trap can then be placed just 
below the atomic trap. 

IV. THE ARPA YIELD DEPENDENCE ON THE PHASE-SPACE DENSITY 
AND COMPARISON TO MAGNETO-ASSOCIATION. 

Following Refs. [H [3], we now present a detailed calculation of ARPA efficiency in a 
thermal ensemble. In order to estimate the fraction of atoms photoassociated per pulse- 
pair, we multiply P{E), the single collision photoassociation probability at energy E, by the 
number of collisions experienced by a given atom while the pulses are on. This is equivalent 
to averaging over all possible values of to as performed above. 

The number of collisions during the pulses is calculated as follows: at a given energy E, 
the velocity of a given atom is f = (2E/m)^ and the distance traversed by it during a pulse of 
Tiaser duratiou is VTiaser- The cross-sectiou for collision is nb"^ where b is the impact parameter. 
For s-wave collisions, the semiclassical estimate is b = h/2p = h/ {2\/2/mE). Hence, the 
number of collisions experienced by the atom during the two pulses is n = Nnb'^vTiaser/V 
where N is the number of atoms in the trap, and V is its volume. Putting all this together 
we obtain that the fraction of atoms photoassociated per pulse-pair is 



Estimating the photoassociation yield for the case of fiat continuum, we can set P{E) ~ 1 
[3] , and assume that all collisions occur at the temperature of the relative motion Trei = 2T 
|24j . Further, when optimizing the yield for an atomic ensemble, we must choose Tiaser — 
2'Kh/kTrei because the bandwidths of the pump and dump pulses should match the energy 
spread in the ensemble. Thus we obtain 



P{E)7lNTiaser 

Wm^l\2E)h 



(37) 



(38) 



V{2mkTre.{)'^l'^ ■ 
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As pointed out above, a narrow resonance can enhance this fraction by a factor of 0.56, i.e., 
for a narrow resonance f"-°-^^'^'^(T^ = 1.56/(T), with /(T) given by Eq. (38). 

We now consider the yield of magneto-association. In this process a time- varying external 
magnetic field is "swept" in magnitude, thereby moving Eg, the position of the Feshbach 
resonance of interest, to lie below Eth, the onset of the continuum. In this way the Feshbach 
resonance is stabilized to become a "Feshbach molecule" [T214T6] . The magneto-association 
is then followed by a traditional 3-bound states STIRAP [THl US] • 

The efficiency of this scheme is limited by the yield of the first step. In this step two 
atoms may form a molecule if prior to sweeping the magnetic field they are within, ap- 
proximately, Qassociation = (^/2)'^ volumc of phase space from each other [23]. Therefore at 
low to moderate phase space densities, the probability for a given atom to participate in a 
magneto-association process is 

f{T^ — N O association/^ whole (2,171 kT ^3/2 C^'^) 

where Qwhoie = V x {2mkTreiY^'^ is the single-particle phase space volume at temperature 
T, and V is the trap volume. 

Comparing Eqs. (38) and (39) we see that the ARPA yield scales with temperature in 
exactly the same fashion as the magneto-association yield. However, in absolute numbers our 
estimates are that at low to moderate phase space densities the ARPA yield is 1.56 x 87r^ ~ 
120 times higher than the magneto-association yield. These findings strongly suggest that 
an experimental investigation of ARPA at sub-/iK temperatures is warranted. 



V. CONCLUSIONS 



In this paper we have shown that Adiabatic Raman Photoassociation of ultracold atoms 
proceeding via collisional resonances is an efficient way of producing ultracold diatomic 
molecules in deeply bound states. We have done that by replacing the resonance-dominated 
molecular continuum by a discrete set of "effective modes" acting like a set of resonances. 
Though when the scattering resonance width is narrow it covers a smaller region in phase 
space (relative to the case of wide resonances), resulting in a drop of the single collision 
transfer efficiency, this drop is amply compensated for by the (as much as an order of 
magnitude) longer durations at which the photoassociating pulses can effectively act. In 
this way each pair of (pump and dump) laser pulses can act on more colliding atoms. The 
overall effect is that the narrow-resonances molecular production yield can be as much as 
~ 56% higher than the wide resonances yield. For atomic temperatures in the /zK range, 
we find that the optimal conditions are attained for resonances whose widths are about 
8 % of the ensemble temperature. We have also shown that the efficiency of the ARPA 
scheme compares favourably with the efficiency of magneto-association, with the yields of 
both schemes scaling with temperature in exactly the same manner. We have demonstrated 
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that the ARPA process is a projective quantum measurement by the pulses of the initial 
continuum wave packet. This feature is a result of the single collision transfer efficiency 
being proportional to the degree of overlap between a function set by the pulses and the 
initial wave packet. 

Future applications will deal with time-dependent resonances. We envision combining 
ARPA with a dynamical sweep of the Feshbach resonance across the threshold energy range. 
As the sweep will render the resonances narrower, the laser pulses will be made narrower so 
as to transfer the atomic gas into molecules in an optimal piecewise manner. 



Appendix: ARPA as a projective measurement of the initial continuum wave 
function. 

We showed in reference [2] that if the continuum is flat then ARPA implements a projec- 
tive measurement of the initial wave function of two colliding atoms. Basically, the profiles 
of the laser pulses e„(t) define a wave form f(^^P^) that is adiabatically coupled to the 
target state |1). An initial scattering state which overlaps f'^^^^'^\t) well will undergo pop- 
ulation transfer to |1), while a state orthogonal to will not. By controlling the laser 
pulse profiles and implementing ARPA one is essentially measuring the wave function of the 
colliding atoms. 

In this Appendix we extend the treatment to a resonance-dominated continuum, and 
relate the effect of collision resonances with our ability to control f'^^^^'^\t). We start by 
noting that in the adiabatic limit the solution of the equations of motion ( 12p3 ) is of the 
form [H E] 



b^{t) = lit) [ dt'exp i [ £+{t")dt" 

Jo L Jt' 



— ism.6{t) I dt'exp 
Jo 

b2{t) = ism9{t) / dt'exp 
Jo 



£^{t")dt" 
S+{t")dt" 



Sme{t') fsourceit') 
) fsourceit ) 
Sme{t') fsourceit') 



(40) 



+ lit) f dt'exp i [ S4t")dt" 

Jo L Jt' 



) fsourcei^ ) • 



where 



and 



s^it) = \ ^T^fsit) ± y/^TM^yp^rvK^)^} 



tan^(t) = £+it)/Vti2it) . 



(41) 



(42) 



(43) 



The final yield of the ARPA process is defined as the probability Pi = |6i(t — )■ oo)\'^. 
Using Eqs. (42) and (43), we see that (t — )■ oo) = 0, and the excited bound state amplitude 
h2it — )■ oo) indeed vanishes. Substituting (t — t- oo) = and sin6'(t — > oo) = 1 in Eq. (40) we 
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obtain that 



POO 

hlit Oo) = fARPAit) fsource(t)dt = {f ^IpaI f source) t (44) 

^0 



where 

r rt 

/arpa(^) = -*exp 





£_(t') dt' 



it). (45) 



Thus the photoassociation amphtude bi{t — )■ oo) is given as the projection of the source 
function f source onto the specific wave form /^Rp^ whose shape is controlled by the amplitudes 
and the phases of ^nit) and ^lp{t). Wave packets that are orthogonal to /arpa do not 
photoassociate in the ARPA process, while the ones that project well onto /arpa do. By 
tailoring the amplitudes and phases of the laser pulses, one can choose which continuum 
waveform is transferred into the target state 

In the main part of the paper we have shown that resonance dominated ARPA is most 
efficient when the resonance is narrow and the pump pulse has a narrow bandwidth relative 
to the initial ensemble temperature. These arguments allow us to replace in (|9| E^h by —oo 
and obtain, via the convolution theorem, 

/oo 
WiT)foit-T)dT , (46) 
-oo 

1 r°° 

W{t) = — ii2AE)e'^^'dE, (47) 



where 



and 



/oo 
6^,(0) e*^-*dE. (48) 
-oo 

/o(t) is the phase-space envelope [21] of the initial wave packet of continuum states. Semi- 
classically it corresponds to the incoming wave function as a function of time t, measured 
at the turning point of a classical trajectory of energy Eq. Positive values of r - the time 



variable in the convolution integral of Eq. (46) - correspond to an outgoing motion, and 
negative values - to an incoming motion [2T1 . 

The function W{t)^ the "FC window", describes the residence time of the system in the 
"FC region," the spatial region for which the FC factors are substantial. The temporal 
width Atvk at which W{t) is substantial, corresponds to Ai?, the spatial extension of the 
FC window. In the flat continuum case, W{t) = /i5(r), and the time of residence in the FC 
region is zero. In contrast, a narrow resonance can cause the system to be greatly delayed 
in the FC window. In that case, an incoming continuum wave packet does not leave the 
Franck-Condon region right after entering it, but rather dwells there for the time given by 
the width Atw- 



We now consider the role of the continuum structure. Combining Eqs.(46) and (47) with 



(44), and introducing t' = t — r we obtain 

POO 

hit ^ oo) = / fZLit')f^'\t') dt' ^ {fZL\f^% , (49) 
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where 



/oo 
jfUit' + T)W{T)dr. 



(50) 



Thus the shape of the wave form scooped from continuum by a laser pulse pair is equally 
defined by the laser pulses and by the continuum structure encoded in the window function 
W{t). 

Expanding the bound continuum transition spectrum into effective modes (Eq.([l4))) we 
have 



W{t) 



exp 



i{E2 - a;2 - -Es) - Y I r 



(51) 



where 6'(r) is the Heaviside function. Therefore 



(52) 



Equations ( 44p2 ) present the main result of this section. They show that, similar to the case 
of a flat continuum, the profiles of the pump and dump pulses define the shape of coherent 
wave forms which can be transferred from the continuum into the target state. However the 
dwelling of the wave function due to the resonances decreases the ability to control these 
wave forms. If due to the resonances the dwell time Ar^y exceeds the durations of the laser 



pulses, then by Eq.(50) we know that in addition to photoassociating atoms that arrive at 



the PC at to, there is a non-negligible probability to photoassociate atoms which get there 
before or after to- 
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